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We focus on the HR of massive vector (spin-1) particles tunneling from Schwarzschild BH ex¬ 
pressed in the Kruskal-Szekeres (KS) and dynamic Lemaitre (DL) coordinates. Using the Proca 
equation (PE) together with the Hamilton-Jacobi (HJ) and the WKB methods, we show that the 
tunneling rate, and its consequence Hawking temperature are well recovered by the quantum tun¬ 
neling of the massive vector particles. 


I. INTRODUCTION 


One prediction of the theory of general relativity (GR) devised by Einstein involves BHs. In principle, a BH is 
classically defined by an area of space called the ’’event horizon”, where everything is swallowed. Beyond the event 
horizon, matter and light flow freely, but as soon as the horizon’s intangible boundary is crossed, matter and light 
become trapped. So, a BH is an invisible object (i.e., black), at least classically. 

Stephen Hawking’s prediction that a BH might not be completely black is unarguably one of the important conse¬ 
quences of the quantum mechanics, when integrated with GR [iHsj. In particular. Hawking proved that a semiclassical 
BH possesses a characteristic temperature of a thermally distributed radiation spectrum, which is the so-called Hawk¬ 
ing radiation (HR) [I|. Today, in the literature there exists several derivations of the HR, which are proposed to 
strengthen this staggering theory (see, for example, IMl). Among those methods, the quantum tunneling method 
(QTM) of Angheben [T^ and Padmanabhan [3,13 (with their collaborators) has garnered much attention (see 
and references therein). QTM employes the complex path i nteg ral analysis of Kerner and Mann Bin in the HJ 
formalism, which takes account of the WKB approximation [3|- According to the QTM, a wave propagator that is 
proportional to exp (■|5'o + S'! -I- 0{h)) is applied to the wave equation of the tunneling particle under question. Here, 
each S denotes the classical action of the trajectory of the particles coming out/in from the horizon. 

In particle physics, a vector boson is a boson with the spin-I. In particular, the massive vector bosons [3 i-6-, 
and Z particles (force carriers of the weak interaction) play a prominent role in the confirmed Higgs Boson [3 - 
Nowadays, the detection of a massive photon, which is the so-called Darkliqht has become very popular in 

the experimental physics since it is envisaged to explain the dark matter [25|. Furthermore, in theoretical physics, 
HR of the massive vector particles in stationary BHs have also attracted much attention (see, for example, (26l - K^ ). 
However, the number of studies regarding the HR of the spin-I particles from the non-stationary regular metrics is 
very limited [3 - s-iid hence those regular spacetimes deserve more research. Such an extension is one of the goals of 
the present paper. For this purpose, we consider the PE [3113,1111 in the KS [3,113 nnd DL M\ coordinates. Next 
we apply the QTM to the PEs, and obtain a set of differential equations for each coordinate system. Those equations 
enable us to get a coefficient matrix. After setting the determinant of the coefficient matrix to zero, we get the action 
So, which is the leading order in h. Then, we show how one can compute the tunneling rate of the vector particles in 
the non-stationary metrics, and recover the standard Hawking temperature of the Schwarzschild BH. 

The paper is organized as follows: In Sec. 2, we first give a brief introduction about the Schwarzschild spacetime in 
KS coordinates. Then, a detailed calculation of quantum tunneling of spin-1 particles near the KS horizon is provided. 
Section 3 is devoted to the computation of the HR of the Schwarzschild BH from the tunneling of the massive vector 
particles in the DL coordinates. The PE with minimum length effect and conclusions are presented in Sec. 4. 
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II. QUANTUM TUNNELING OF MASSIVE VECTOR PARTICLES FROM SCHWARZSCHILD BH IN 

KS COORDINATES 

The well-known Schwarzschild solution is in general described by the coordinates t and r as follows 

ds^ = gttdf + grrdr'^ + geedd'^ + g^^dip"^, ( 1 ) 

where gu = (1 — 2M/r), g^r = —(1 — 2M/r)~^, ggg = —r^, and g^^ = —r^ sin^ 6. Herein, there is an event horizon at 
r — 2M, so that grr blows up. On the other side when r < 2M, the gtt and grr exchange their signatures, however the 
signatures of ggg and g^p^p are not affected. Hence, r becomes “timelike” and t becomes “spacelike” inside the event 
horizon. One could clear up this “coordinate singularity” problem by introducing the KS coordinates [i^, : 


ds^ = Ai^—dr'^ -I- dB?) + r^{d9'^ + B^dip^), 


( 2 ) 


where B = sin(0) and the metric function A is given by 


A = 


32M^ 

r 


(3) 


Metric ([2]) covers the entire spacetime manifold of the maximally extended Schwarzschild solution, and it is well- 
behaved everywhere outside the physical singularity (r = 0). The event horizon in the KS coordinates corresponds to 
T = ± i?, and the curvature singularity is located at = 1. Furthermore, in this coordinate system the Killing 

vector becomes 


e = 


—,—. 0.0 

AM’ AM' ’ 


The particle energy of a test particle is given by (in terms of the action S) [^, \A 




On the other hand, for a curved spacetime, the PE is governed by 


(4) 

(5) 


where 'I'y 
and 


1 

7^ 


_ 'ryi 


( 6 ) 


and m represent the spinor fields and mass of the spin-1 particle, respectively. 


Using metric m in Eq. (jG]) , we obtain the following set of differential equations: 


(7) 


- (HArB) ^ -H h^r'^B'^{-diiR^o + d^R^i) - h^B'^A{deeAiQ - dgrA> 2 ) 

+ h?BA {dgB) {drA!2 - dgAio) - h^A{dppAio - = 0, (8) 


(HArB) ^ [m^Ar^B^A/i — h^r^B^{dTRA>o — OttA'i) — h^B^A{deeA'i — 9e_R4'2) 

+ h^BA {dgB) {dRA >2 - dgA>,) - h^A{d^^Afi - = 0, (9) 


(ji'/Ar'^B^ [irL^Ar'^B‘^Ai2 — h^r‘^B‘^{dgrA/o - drTA>2) - B^ {dguA/i - diiRA>2) 

- h'^A{d^^A '2 - dg^A/ 3 )] = 0 , ( 10 ) 
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(hy/Ar^Bi^ ^ [m^Ar^B-^s - fi^r^B{d^r'i'o - drr'^ 3 ) + h^r^Bid^R-^i - dRR'Ss) 

+ h^AB{de^^2 - Seed's) + h^A {deB) {ds'fs - = 0. (11) 

Applying the WKB approximation (s^ : 


'I'jy = exp 

and taking the lowest order of h, Eqs. (8-11) become 


-5'o(r, R,9,ip) + Si{T,R,e,ip) +0{h) 
n 


( 12 ) 


{dRSof - A (d^Sof - m^Ar^B^ - AB^ (dgSof) cq 

+ C 3 A {dr,pSo) + cir^B^ {dRrSo) + C 2 AB‘^ {dreSo) = 0, (13) 


r'^B^ {dRrSo) Co + Cl (-r'^B^ {OrS^f + A (d^Sof + m'^Ar^B^ + AB^ {dgSof) - C3A {Or^So) 

-C2AB^{dReSo) = 0, (14) 


r^B^ [(a.eS'o) CO - Cl (^fleS-o)] + C2{AB^ {d^Sof 


- B 


4^2 


{drSoT - (dRSor - } - AB^cs (d^gSo) = 0 , (15) 


r'^B [(5r^5'o) Co - ci (S^^S'o)] - C 2 AB (d^gSo) 


+ Br^ 


2 ^ 


m^A - {drSoY + (dRSor - ^ (deSo) 


C 3 = 0 . (16) 


Now, one can obtain a matrix equation Z (cq, ci, C 2 , C 3 ) = 0 


transposed vector, and Z represents a 4 x 4 matrix) with the following non-zero elements: 


(the superscript T means the transition to the 


Zii — 

Z 12 = 
^13 = 
Zi4 = 

Z 22 = 
Z 23 = 
Z 24 = 
Z 33 = 

Z34 = 
Z44 = 


-r^B^ (dRSof - A (d^Sof - m^Ar^B^ - AB^ ( 9 e^o)"J , 

Z21 = YB^ (dRrSo) , 

AB^ (drsSo) , Z31 = r^B^ (drffSo) , 

A (dr^So) , Z41 = r^B (dr^So), 


-r^B^ (drSo) + A (d^SoY + vYAYb'^ + ^dgSoY 


-AB^ (dRoSo), Z 32 = -r^B^ (dRgSo) , 

—A (dR^So ), Z42 = — Yb (dR^So ), 

AB^ {d^Sof - B^r^ [(a.^o)" - (a^^o)' - rn^A 

—AB^ (dipgSo ), Z43 = —AB (ai^gS'o), 


BY 


rYA - (drSoY + (dRSof - ^ (dgSo) 


Let us consider the following HJ solution for the action: 


(17) 


'S'o = Q(t, R) + k (9) + j(f, 


(18) 
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where j denotes the angular momentum of the massive vector particle. Thus, the determinant of Z-matrix yields 


detZ = -m^r^AB'-^{A {dgkf + f 




rn^A - [drQ f + {duQ^ } 


(19) 


The nontrivial solution for dnQ is obtained by the condition of ”detZ = 0” [1^. Hence, after substituting drQ = 
~ [recall Eq. (0] into Eq. (19), we obtain 


= 


AEMtvB ± R^lQE'^M^r^B^ - A{R^ - r^) | {dgkf + B^ + j2| 


(i?2 - r2)Hr 


( 20 ) 


by 


where + (—) corresponds to the outgoing (incoming) massive vector particles. The definite integration of Q is given 


Q = J {9rQ) dR + {drQ) dr. 


( 21 ) 


Using the identity drQ = — ^d^Q, once again, Eq. (1^ can be rewritten as 


Q=l [^-^d{R^-rf-^-^ 
^ 2 J R ^ ^ R 


dr. 


( 22 ) 


It is obvious that the second term is real in Eq. (1^^ . However, after inserting Eq. (uni) into Eq. (1^ . we see that 
the imaginary contribution to the action comes only from the first term since it has pole at the horizon. Thus, the 
complex path integration method [H, [l^ for the pole located at the horizon {R = r) yields 


ImQ- 


horizon 


= 0 , 


(23) 


IvtlQ-\-\horizon — 4l7T AIE, 

Therefore, the probabilities of the ingoing/out going massive vector particles become 

p — p—-f[^'mQ-\horizo-n. — 1 

^ absorption — ^ 


(24) 


(25) 


r 


emission 


_^-^ImQ+\horizon — g-S-n-MB 


(26) 


It is worth noting that the above results are in full agreement with the semiclassical QTM [13, which expects a 
100% chance for the ingoing particles to enter the BH, i.e., Tabsorption = 1; and thereupon computes the probability 
of the outgoing (tunneling) particles, Temission- 
The tunneling rate is then computed by 


r = 


r 


emission 


r absorptii 


g —SttME 


(27) 


Now, recalling the Boltzmann factor (see for example 113), r = e = e where /3 is the inverse temperature 

we can recover the original Hawking temperature of Schwarzchild BH: 


T = Th 


fin) 

47r 


1 

8ttM 


(28) 
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III. QUANTUM TUNNELING OF MASSIVE VECTOR PARTICLES FROM SCHWARZCHILD BH IN 

DL COORDINATES 

In 1933, Georges Lemaitre found a coordinate system ^f, R, 9, that removes the coordinate singularity at 
the Schwarzchild BH is given by 


~F 


ds^ = -dr^ + -+ AM^F\d9^ + B^dp^) 


(29) 


where 


F = 




(30) 


The event horizon in the DL coordinates corresponds to F = 1 or i? = ^ + f. 
reads 


Furthermore, the Killing vector 


e= [1,1,0,0], (31) 

and it leads to the following particle energy (i^ : 

E = -ed^S = - {dr- + Or- ) S. (32) 

In this coordinate system, PEs ([S]) with the ansatz (TT^ are given by 

Co (- {d^Sof - AF^Ad^B^ (9^5'o)^ - Am^F'^M'^B'^ - B^ {deSof) + AciF^M^B^ {dpifSo) 

+ C2B^ (dfeSo) + C 3 {dfipSo) = 0, (33) 

AF^M^B'^co {df^~So) + ci (Am^F^M^B^ - AF^M^B^ {dfSof + B^ {dgSof + (a^Fo)") 

— C 2 B (dj^gSo) — C 3 = 0, (34) 


AF'^Ad^B^co {dfeSo) - AciF'^M^B^ci {dj^gSo) 

+ C 2 [Atu^F^M^B^ + [d^Sof - AF^Ad^B^ {dfSof + AF^M^B^ [dRSof] - C 3 (a^gFo) = 0, (35) 


AF^Ad^Bco (df^So) - AF^Ad^Bci (a^^ 5 o) - C2B (d^eSo) 


Am^F^Ad^B + B {deSoY - AF^Al^B {dfSoY + AF^Ad^B {dp^SY C 3 = 0. (36) 


Now, one can read the non-zero elements of the coefficient matrix H (co, ci, C 2 , C 3 ) = 0 (H is another 4x4 matrix) 

as follows 
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Hii 

Hi2 

Hi3 

Hi 4 

^22 

^23 

^24 

^33 

H 34 

H44 


- {dy^Sof - AF^M^B^ {diiSof - Am'^F'^M^B^ - B^ (ae 5 o)^J , 
AF^M^B^ {dji-So), H21 = AF^M^B^ (dj^-So) , 

B^ (dreSo) , H31 = AF^M^B^ (d^eSo) , 

(df^So), H41 = AF^M^B (df^So) , 

Am'^F^M'^B'^ - AF'^M^B^ {d^Sof + {dgSof + ( 5 ^S'o)^j , 

-B^ {dj^gSo), H32 = -AF^M^B^ [dp,,So) , 

- {dR^So) , H42 = -AF^M^B (df^^So) , 

[Am^F'^M^B^ + {d^Sof - AF^M^B^ (dfSo)^ + AF^M^B^ (9flS'o)^], 

— (d^eSo), H43 = —B (dipgSo), 

Am^F^M'^B + B (dgSo)^ - AF^M^B (dfSo)^ + AF^M^B (dj^So)^ . 


Inserting the following ansatz for Sq 


So = Q{f, R) + k{e) + jip, 

into Eq. (1571) . and subsequently using the energy condition (|57|) . namely: 


(37) 

(38) 


we get solutions for df^Q from detH = 0: 


dfQ — —{E F dj^Q), 


(39) 


detH = — 


M^F^Bm^ 

1024 


B^ (dgkf + f + AM^F'^B'^ 


m + F 




= 0 , 


as follows 


EMBF ± ^jE^M'^E^B'^ - {F - 1) m^F'^M'^B^ + \B‘^ {dgk)^ + ^ 


(F - 1) FBM 

Using the energy expression (IMl) in the definite integration of Q : 


(40) 


(41) 


Q = J df^QdRF dfQdf, (42) 

we obtain 

Q — J dpQd (^R — — E j dr, 

= 2M J ^^dE - E 

where dF = -^VFd{R — f) [see Eq. ([30|) ]. It is clear that the second integral of (|43l) results in real values, which 
means that it does not give any contribution to the imaginary part of the action. However, after substituting Eq. 
(141) into Eq. (|43)) . one can see that the first integral of Eq. (1431) has a pole at the horizon (F = 1), and evaluating it 
around the pole yields 



lTnQ+\horizon — AtT M E ^ 


(44) 


and trivially 


ImQ -1 

horizon — 0- (45) 

The above results are fully consistent with Eqs. ([23ll and (l24ll . Consequently, they admit the same tunneling rate 
given in Eq. (1271) . In short, using the quantum tunneling of the massive vector particles in the DL coordinate system, 
we have managed to rederive the Hawking temperature [Th = -^^m) of the Schwarzchild BH. 
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IV. CONCLUSION 


In this paper, we have used the PE (|6l) in order to compute the HR of the massive vector particles tunneling from 
the Schwarzchild BH given in two different (KS and DL) regular dynamic coordinate systems. In addition to the 
HJ and the WKB approximation methods, particle energy definitions played crucial role in our computations. The 
original Hawking temperature of the Schwarzschild BH is impeccably obtained in the both coordinate systems. Thus, 
we have shown that HR is independent of the selected coordinate system. The latter remark was also highlighted in 

sa. 


Finally, we plan to extend our present study to the HR of the massive spin-1 particles that experience the minimum 
length effect , which is governed by the GUP (generalized uncertainty principle) . Such a study requires 

the following modification in the PE 


1 

7 ^ 







m 

1 ^ 


= 0 , 


where 


(46) 




-d. 


' 0 


d 


(47) 


in which Ip denotes the Planck length. This problem may reveal more information compared to the present case. This 
is going to be our next study in the near future. 
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